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Abstract 


Micrometeorological  field  data  were  used  in  an  effort  to  evaluate  an  atmo¬ 
spheric  microclimate  model,  specifically,  to  determine  if  we  can  forecast 
one-dimensional  profiles  of  the  atmospheric  boundary  layer  sufficiently 
well  that  the  corresponding  outdoor  acoustic  field  can  be  reliably  calculated . 

A  soil-plant-atmosphere  model  was  exercised  with  two  sets  of  experimen¬ 
tal  data.  The  model  was  used  to  calculate  the  surface  energy  budget  and  to 
derive  meteorological  profiles  for  the  overl5nng  boundary  layer;  both  of 
these  were  compared  to  field  data.  These  modeled  results  were  also  used  to 
derive  outdoor  sound-speed  profiles  and  were  used  as  input  to  a  short-range 
acoustic  propagation  numerical  code. 
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1.  Introduction 

In  the  Glossary  of  Meteorology  (Huschke,  1959),  "microclimate"  is  defined  as 
"the  fine  climate  structure  of  the  air  space  which  extends  from  the  very 
surface  of  the  earth  to  a  height  where  the  effects  of  the  immediate  character 
of  the  imderlying  surface  no  longer  can  be  distinguished  from  the  general 
local  climate  (mesoclimate  or  macroclimate)."  In  modeling  microclimate, 
one  simulates  the  characteristics  of  the  atmospheric  surface  and  boundary 
layers.  Its  time-varying  behaviors  are  related,  albeit  nonlinearly,  to  radia¬ 
tive  heating  and  cooling,  changes  in  water  content  of  both  the  air  and  soil, 
terrain,  land  use,  and  ground  cover.  Numerical  computer  models  attempt 
to  simulate  the  microphysical  processes  of  microclimate  for  a  wide  range 
of  applications.  In  defense  technology,  microclimate  data  can  be  applied 
operationally,  as  well  as  in  support  of  planning,  environmental  assess¬ 
ment,  and  research. 

Several  turbulence-diffusion,  atmospheric  boundary-layer  models  are  pre¬ 
sented  in  the  literature,  including  those  of  Avissar  et  al  (1986),  Burk  (1977, 
1980),  Blackadar  (1978, 1979),  Cionco  (1965, 1985),  Deardorff  (1972, 1978), 
Mellor  and  Yamada  (1974),  Naot  and  Mahrer  (1989),  Pielke  and  Mahrer 
(1975),  Yamada  and  Mellor  (1975),  and  Zhang  and  Anthes  (1982).  With 
some  proven  measure  of  accuracy,  they  attempt  to  determine  the  temporal 
changes  of  specific  meteorological  quantities  for  different  environments 
and  sets  of  initial  conditions.  Many  of  these  models  represent  sets  of  equa¬ 
tions  for  two-  and  three-dimensional  mesoscale  (20  to  200  km  domain)  pre¬ 
dictions  over  complex  terrain. 

This  report  documents  an  evaluation  effort  in  which  I  studied  a  one¬ 
dimensional  (ID)  version  of  the  plant-canopy  microclimate  model,  as 
given  by  Avissar  and  Mahrer  (1982,  1988),  using  two  dissimilar  sets  of 
micrometeorological  field  data  (Clarke  et  al,  1971;  Stenmark  and  Drury, 
1970).  I  made  model  calculations  of  the  surface  energy  budget  and 
compared  them  to  field  observations  for  both  data  sets.  I  also  computed 
meteorological  profiles  for  the  overlying  boundary  layer,  comparing  these 
also  to  the  observed  field  data.  I  used  these  modeled  results  to  derive  out¬ 
door  sound-speed  profiles.  The  model  output  data  were  also  used  as  input 
to  the  acoustic  propagation  numerical  code  WSCAFFIP  (Noble,  1996; 
Noble  and  Marlin,  1995)  for  calculating  levels  of  short-range  acoustic 
attenuation.  My  goal  was  to  determine  if  it  was  possible  to  forecast  ID 
profiles  of  the  atmospheric  boundary  layer  sufficiently  well  that  the 
corresponding  outdoor  acoustic  field  could  be  reliably  calculated.  Time- 
dependent  forecasts  of  microclimate  processes  are  rarely  used  in  atmo¬ 
spheric  acoustics  (Noble,  1991, 1992). 

Following  an  overview  of  the  soil-plant-atmosphere  model  (sect.  2),  I 
describe  the  two  evaluation  data  sets  (sect.  3).  Model  results  are  discussed 
in  section  4. 1  then  present  the  derived  sound-speed  profiles  and  results 
from  the  acoustic  propagation  numerical  code  WSCAFFIP  (sect.  5).  Finally, 
conclusions  are  drawn  (sect.  6)  and  recommendations  related  to  this  area 
of  study  are  given  (sect.  7). 
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2.  Model  Description 

The  ID  turbulence-diffusion  model  of  the  plant-canopy  microclimate 
reported  by  Avissar  and  Mahrer  (1988)  is  based  on  their  earlier  efforts 
(Avissar  and  Mahrer,  1982),  as  well  as  the  earlier  works  of  Deardorff  (1978) 
and  Pielke  and  Mahrer  (1975).  It  is  a  simulation  model  that  generates  a 
fairly  complex  set  of  calculations  at  the  soil  and  plant  level,  so  that  the 
surface  energy  budget  (e.g.,  the  transference  of  heat  and  moisture)  can  be 
made  integral  to  processes  within  the  first  1  to  3  km  of  the  overlying 
atmosphere.  Although  the  model  was  originally  developed  for  applica¬ 
tions  in  irrigation  management  by  engineering  estimates  of  actual  and 
potential  evapotranspiration  (Avissar  et  al,  1986;  Kordova  et  al,  1994),  I  re¬ 
vised  it  to  yield  time-dependent  calculations  of  the  atmospheric  acoustic 
field  that  result  as  the  derivative  of  modeled  wind  speed,  wind  direction, 
temperature,  and  humidity  profile  estimates. 

As  initial  input,  the  model  requires  the  numerical  equivalent  of  the  day  of 
the  year  and  the  time  of  day  (relative  to  the  Greenwich  meridian),  site  loca¬ 
tion  data  (latitude  and  longitude),  fraction  of  sky  cloudiness,  and  ground 
cover  data  (i.e.,  canopy  height,  leaf  area  index,  surface  reflectivity  (albedo), 
thermal  emissivity,  and  aerodynamic  roughness  length).  The  model  com¬ 
putes  soil  properties  (e.g.,  hydraulic  and  thermal  conductivity  and  soil 
specific  heat  capacity)  from  inputs  of  soil  water  content,  porosity,  texture 
(i.e.,  proportion  of  sand,  clay,  and  organic  matter),  and  bulk  density.  Other 
model  constants  include  subsoil  properties,  such  as  plant  root  density  and 
distribution.  Also,  initially,  the  model  requires  profile  data  for  air  tempera¬ 
ture  and  specific  humidity  (atmospheric  water  vapor  content)  from  the 
ground  level  to  the  model  top.  Lastly,  the  screen  level  (2  m  above  ground 
level)  values  of  wind  speed,  wind  direction,  temperature,  and  humidity 
are  updated  at  each  hour  of  the  model  run.* 

The  model's  estimates  of  ID  profile  structure  are  computed  from  a  simpli¬ 
fied  set  of  kinematic  and  thermodynamic  equations.  The  transference  of 
momentum  (i.e.,  effects  of  surface  roughness  and  wind  shear)  and  the  ex¬ 
changes  of  heat  and  moisture  from  one  level  to  the  next  are  calculated  in 
terms  of  eddy  coefficients  and  vertical  gradients  that  represent  the  flux. 
The  model  turbulence-diffusion  equations  for  momentum  can  be  ex¬ 
pressed  as 


du  r  r 


jj)  .  (2) 

for  the  u-  and  Z7-’Components  of  wind  speed,  where/ denotes  the  Coriolis 
parameter  (the  deflecting  force  due  to  the  earth's  rotation  acting  upon 

*  Although  this  model  requires  a  tremendous  amount  of  initial  input  as  well  as  supporting  data  for  updating  calcula¬ 
tions  every  hour,  it  may  still  be  possible  to  apply  it  generally  to  locations  worldwide.  One  could  use  bulk  estimates  of 
the  soil  thermal  properties  and  plant  or  ground  cover  data,  and  hourly  meteorological  inputs  (such  as  wind  speed, 
temperature,  or  cloud  cover)  could  be  extracted  from  3D  mesoscale  model  forecasts  run  independently  from  the  micro¬ 
climate  code. 


^  dz\  dz 


(1) 
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moving  air),  the  subscript  g  refers  to  the  geostrophic  wind  (such  that  the 
second  terms  on  the  right  represent  the  pressure  gradient  force),  and 
denotes  the  exchange  coefficient  for  the  diffusion  of  momentum.  In  the 
surface  layer  Km  is  calculated  as  Km  =  ku^f  where  k  is  Karman's  con¬ 
stant  (a  dimensionless  constant  {k  =  0.4)  of  atmospheric  fluid  dynamics),  z 
is  height  (in  meters),  m*  is  the  friction  velocity  (in  meters  per  second)  that 
quantifies  drag  or  surface  stress,  and  ^  is  the  nondimensional  wind  shear, 
which  accounts  for  stability,  turbulence,  and  near-ground  wind  speed  pro¬ 
file  structure  (Monin  and  Obukhov,  1954;  Businger  et  al,  1971;  Dyer,  1974). 

In  the  convective  (heated),  daytime  boundary  layer,  the  model  derives  the 
mixing  coefficients  as  a  function  of  height,  as  suggested  by  O'Brien  (1970). 
The  function  can  be  expressed  as 


k,=k,M 


'  (Zj-zY 

i(2,-2s)' 


8K, 


dz 


(3) 


where  z  is  height  above  ground  level  and  the  subscripts  s  and  i  refer  to  the 
heights  of  the  top  of  the  surface  layer  and  top  of  boundary  layer,  respec¬ 
tively.  The  profile  function  K^  is  applied  similarly  for  the  heat  and  moisture 
coefficients  where  for  heat  exchange  K^  =  kd^zl  and  for  moisture 
exchange  K^  =  for  z  <  Zj,  where  and  are  scaling  variables  that 

quantify  the  near-ground-level  gradients  (Rachele  et  al,  1995,  1996a)  of 
temperature  and  specific  humidity,  respectively.  The  model  defines  Zj  = 
0.04  z,  and  sets  Kz,  =  1-0  for  z  >  z,. 

Finally,  the  growth  and  decay  of  the  planetary  boundary  layer  is  calculated 
from  an  expression  derived  by  Deardorff  (1974),  which  can  be  stated  as 


l.sjiv!  +  l.lul-33utfz^ 


(4) 


where  iv^-  is  the  vertical  velocity  at  z,  (assumed  to  be  negligible  or  equal  to 
zero,  since  time-dependent  calculations  of  vertical  velocities  are  not  explic¬ 
itly  derived),  Zg  is  the  roughness  length  that  characterizes  the  frictional  ef¬ 
fects  of  the  land  morphology  and  the  type  of  ground  cover,  d6P/ dz  is  the 
vertical  gradient  of  potential  temperature  (pressure  scaled  temperature)  in 
the  stable  air  immediately  above  the  boundary-layer  top,  and  finally  the 
vertical  velocity  scaling  variable,  an  implicit  calculation  of  buoyancy,  is 
defined  as 

w.  =  [{-g/e)u,6^z]^'^  (5) 


where  g  is  acceleration  due  to  gravity. 

The  set  of  equations  for  calculations  through  the  soil  layer,  surface  layer, 
and  plant  canopy  are  also  given  by  Avissar  and  Mahrer  (1988).  However, 
in  this  work  I  made  substitutions  based  on  previous  work  of  Rachele  and 
Timick  (1994)  for  the  equations  used  to  determine  solar  angles  (Woolf, 
1968),  insolation  (Meyers  and  Dale,  1983),  and  clouds  (Haurwitz,  1945). 
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The  surface  energy  budget  is  applied  to  three  layers:  the  soil  (including 
several  thin  sublayers);  the  vegetative  canopy  and  the  air  throughout  the 
canopy;  and  the  thin  (1  to  4  m)  layer  of  air  that  extends  above  the  canopy 
top  (Avissar  et  al,  1986).  Air  temperature  and  water  vapor  content  for  each 
layer  are  derived  as  a  result  of  energy  and  mass  transfers  within  the  sys¬ 
tem  (including  plant  transpiration).  TT^e  model  also  includes  heat  conduc¬ 
tion  and  moisture  diffusion  equations  for  the  soil  layer,  since  soil  wetness 
and  surface  temperature  are  possibly  the  two  most  significant  factors 
determining  the  partitioning  of  the  total  available  energy  to  the  soil, 
evaporative,  and  sensible  heat  fluxes.  The  total  available  energy  (i.e.,  the 
net  radiative  flux)  is  approximated  as  a  function  of  transmission,  terrain, 
land  use,  albedo,  canopy  leaf  area  characteristics,  soil  wetness,  and  soil 
thermal  properties.  The  equations  of  the  plant-canopy  energy  budget  are 
given  in  table  1. 

Table  1.  Components  of  atmospheric  surface  layer  and  plant-canopy  radiation  and  energy  budget 

(Avissar  and  Mahrer,  1988). 


Component  Equation  Etefinition  of  terms 


Plant-canopy  energy 
budget  (Carson,  1987) 

Rj^  =  net  radiative  flux 
=  sensible  heat  flux 
=  evaporation  heat  flux 

Sq  =  soil  heat  flux 

Short-wave 
radiation 
(insolation),  Rg 
(Meyers  and  Dale, 
1983) 

Rs  =  kTRTGTwTA  COS  (Z) 

model 

top 

4  =  solar  constant  (=1367  W/m“^) 
Ti  =  transmission  coefficients 
for  Rayleigh  scattering  (R), 
absorption  by  permanent 
gases  (G)  and  water  vapor 
(W),  and  scattering  due  to 
aerosols  (A) 

Z  =  solar  zenith  angle  (°) 

Sky  long-wave 

a  =  Stefan-Boltzmann  constant 

radiation,  Rii 

2  =  0 

(5.6697  xl0-«Wm'2K-^) 

Tz  =  air  temperature  (K) 

^  =  ^H20  +  ^02  =  total  emissivity 
%20  ~  emissivity  of  water  vapor 
eco2  ~  emissivity  of  carbon  dioxide 

Net  radiation 

Rnv= 

y/f  =  shielding  factor  (full 

of  canopy,  Rj^ 

{l  +  all-Yf{l-t,)]}Rs+¥fe^ 
([1  +  (1  -  y/jp  -  e^)]R  a  +  tJ} 

canopy  =  1;  bare  soil  =  0) 

Oj;  =  albedo  of  plant  surface 
ag  =  albedo  of  soil  surface 
(albedo  =  reflectivity) 
fp  =  transmissivity  of  the  canopy 
£p  =  plant  foliage  emissivity 

Eg  =  soil  emissivity 

Tg  =  air-groimd  temperature  (K) 
Ty  =  canopy  virtual  temperature 
that  accounts  for  effects 
of  moisture  in  air  (K) 
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Table  1.  Components  of  atmospheric  surface  layer  and  plant-canopy  radiation  and  energy  budget 
(Avissar  and  Mahrer,  1988)  (cont'd). 

Component  Equation  Definition  of  terms 


Net  radiation 
at  ground,  R^g 


^NG  =  (1  -  «g)(l  -¥f+¥f  QRs 

+ (1  -  Ra  +  ¥f  Sg£v  ^  Ty 


Rj  =  short-wave  radiation, 
insolation  (W/m^) 


Sensible  heat,  Ch 

Sensible  heat 
between 

canopy  and  air,  Chv 

Sensible  heat 
between 

ground  and  air,  C^g 

Evaporation  flux,  Ey 

Evaporation  flux 
between  canopy 
and  air,  Eyy 

Evaporation  flux 
between  ground 
and  air.  Eye 

Soil  heat  flux,  Sq 


Ch  =  Chv  +  Chg 
Chv=  VfpCpU*^ 

CHG  =  ('^-¥f]pCp^^*S* 
Ey  =  Eyy  +  EyG 

£vg=(i- 


Heat  conduction 


Ps~ 


g  _ 


dt 


_a 

dz 


A 


dz 


Soil  moisture 
diffusion 

(Naot  and  Mahrer, 
1989) 


dt 


dz 


Da 


dKh  , 

+  -^-LpT 


, _ flux  of  canopy 

~  total  flux  of  canopy  +  soil  l+2Aiy/j: 
Ai  =  leaf  area  index 
p  =  density  of  air  (kg/ m^) 

Cp  =  specific  heat  of  air  (J/kgK) 

=  surface  friction  velocity  (m/ s) 

0.  =  surface  scaling  temperature  (K) 


L*  =  heat  of  vaporization  or 
condensation  (J/kg) 

=  surface  scaling  specific  humidity 

(g/g) 


A  =  soil  thermal  conductivity 

(W  m"^  K“^;  derived  in  terms 
of  soil  composition  and  soil 
wetness) 

Cs  =  soil  specific  heat  (J  m'^  K"^) 

Ps  =  soil  density  (kg  m“^) 

A  =  soil  water  content  (%) 

Da  =  isothermal  moisture  diffusivity 
Dt  =  thermal  moisture  diffusivity 
Kfi  =  hydraulic  conductivity  (m/s) 
Lpj  =  soil  water  loss  through 
plant  transpiration 
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3.  Field  Data 


3.1  Hay  1967 

One  of  the  sets  of  experimental  data  used  in  this  study  was  collected  dur¬ 
ing  a  field  campaign  near  the  town  of  Hay,  New  South  Wales,  Australia 
(34.50,  144.93),  which  took  place  during  the  1967  southern  hemisphere 
winter,  from  15  July  to  27  August  (Clarke  et  al,  1971).  This  very  large  data 
set  was  collected  over  an  extensive  flat  area  (3600  km^),  where  the  groimd 
cover  was  predominantly  sparse  low  vegetation  over  dry  soil.  For  the  task 
at  hand  (model  verification  with  regard  to  surface-layer  and  overlying 
boundary-layer  profile  structure),  the  data  were  highly  suitable  in  all  areas 
except  for  plant-canopy  interactions.  Five  types  of  surface-layer  and 
mixed-layer  data  were  collected;  (1)  hourly  16-m  micrometeorological 
tower  measurements  of  wind  speed,  specific  humidity,  and  temperature; 
(2)  hourly  surface  measurements  of  the  net  radiative  balance  and  soil  heat 
flux;  (3)  measurements  of  wind  speed  and  wind  direction  from  0  to  2  km 
(AGL — above  ground  level)  obtained  by  hourly  pibal  (pilot  balloon) 
flights;  (4)  measurements  of  pressure,  temperature,  and  mixing  ratio  (m,  = 
masswater  vapor/n^asSdry  air)/  to  3  height  of  2  km,  taken  by  radiosonde  (instru¬ 
mented  balloon  sonde)  flights  at  3-hr  intervals;  and  (5)  reports  of  frac¬ 
tional,  low,  high,  and  total  cloud  cover  (including  cloud  type  descriptors). 

Of  the  1050  hourly  reports  (profiles)  of  published  observations,  I  used 
those  collected  during  the  27-hr  period  from  2100  LT  (local  time)  on  15  Au¬ 
gust  1967  to  2400  LT  on  16  August  1967.  This  data-collection  period  was 
characterized  by  clear  skies  and  relatively  light  daytime  wind  speeds.  Dry 
soil  promoted  strong  surface-based  convective  heating.  The  soils  for  the 
field  area  were  predominantly  brown  loams  overlying  red-brown  clay  sub¬ 
soils. 

3.2  Davis  1966 

The  second  set  of  experimental  data  used  in  this  study  was  collected  dur¬ 
ing  a  field  study  at  Davis,  California  (38.33, 121.44),  during  the  summer  of 
1966  and  spring  1967  (Stenmark  and  Drury,  1970).  The  field  site  was  a  flat, 
fescue-grass-covered  50,000-m^  area  at  17-m  elevation  (above  mean  sea 
level).  It  was  about  2  km  west  of  the  main  portion  of  the  University  of  Cali¬ 
fornia  at  Davis  campus,  113  km  northeast  of  San  Francisco.  The  data  were 
taken  during  periods  when  fields  surrounding  the  test  area  were  for  the 
most  part  crop  covered  and  well  irrigated,  giving  (in  effect)  horizontal  uni¬ 
formity  to  the  surface  conditions  with  respect  to  temperature  and  mois¬ 
ture.  Advection  effects  due  to  the  site's  proximity  to  the  coast  were  consid¬ 
ered  to  be  negligible  during  the  data-collection  period.  Profiles  of  wind 
speed,  temperature,  and  specific  humidity  were  taken  at  nine  levels  from 
25  to  600  cm.  In  addition,  infrared  radiometers,  ground  heat  flux  plates, 
and  a  large  (3.0-m-diameter)  weighing  lysimeter  were  used  to  record  the 
net  radiation,  the  soil  heat  flux,  and  surface  evaporation  rates,  respectively. 
The  energy  budget  equation  was  used  to  determine  the  sensible  heat  flux 
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Ch  as  the  residual.  The  soils  at  the  Davis  field  site  were  predominantly  silt 
loams. 

In  the  present  study,  I  used  the  30-hr  subset  of  data  taken  from  0100  LT  on 
13  July  1966  to  0600  LT  on  14  July  1966.  This  data  interval  was  character¬ 
ized  by  clear  skies,  damp  soil,  low  wind  speeds  through  the  morning,  and 
moderate  wind  speeds  through  the  afternoon.  Increased  evaporation  rates 
associated  with  the  change  (increase)  in  surface  wind  speeds  were 
observed.  (The  surface  wind  direction  (190°  to  210°)  remained,  for  the  most 
part,  constant.)  The  data  collected  from  this  site  were  very  well  suited 
for  verifying  the  model  with  regard  to  surface-layer  energy  budget 
calculations. 


3.3  Model  Parameters  from  Data  Sets 

Table  2  gives  model  parameters  for  both  the  Davis  1966  and  Hay  1967  data 
sets.  Typical  values  given  for  albedo  and  surface  roughness  for  other  com¬ 
binations  of  soils  and  ground  cover  can  be  retrieved  from  Hansen  (1993a, 
1993b).  Hourly  meteorological  input  for  both  data  sets  is  given  in  tables  3 
and  4. 


Table  2.  Model  ,,  .  ,  , 

.  Variable  name 

parameters  for  Hay 

1Qfi7  and  Davis  1966 

15-16  August  1967 
Hay,  Australia 

13-14  July  1966 
Davis,  CA 

data  sets.  Latitude 

34.50  S 

38.33  N 

Longitude 

144.93  E 

121.44  W 

Albedo 

0.20 

0.27 

Emissivity 

0.98 

0.96 

Roughness  length 

0.13  cm 

1.40  cm 

Ground  cover/canopy  height 

barren  field /O.O  cm 

dense  field/ 10.0  cm 

Leaf  area  index 

0.0 

1.0 

Soil  properties: 

Water  content 

0.08  % 

0.20  % 

Thermal  conductivity 

0.539  W  m-i  K-i 

0.947  Wm-iK-i 

Volumetric  specific  heat 

1.374  Jm-^K-i 

1.887  Jm-^K-i 

Density 

1400.0  kg  m"^ 

=1400.0  kg  m-3 

Soil  texture: 

%  sand 

28.0 

20.0 

%  clay 

70.0 

77.0 

%  organic 

2.0 

3.0 
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Table  3. 

Meteorological  model 
input  (1-hr  average 
data)  for  Hay  15-16 

Local 

time 

Temperature 

(°C) 

Specific  humidity 
(g/kg) 

Wind  speed 
(m/s) 

Wind  direction 
(°N) 

21 

5.4 

4.818 

1.75 

165 

August  1967  data  set 

22 

4.0 

4.266 

1.67 

155 

(winter). 

23 

3.3 

4.295 

1.39 

156 

24 

2.8 

4.190 

1.57 

160 

1 

2.2 

4.058 

1.71 

160 

2 

2.2 

4.058 

2.18 

162 

3 

1.7 

3.957 

2.03 

159 

4 

0.6 

3.693 

1.67 

159 

5 

-0.3 

3.495 

0.77 

158 

6 

-0.3 

3.531 

0.40 

155 

7 

-0.4 

3.506 

0.40 

142 

8 

0.9 

3.774 

1.21 

138 

9 

5.3 

4.893 

2.71 

107 

10 

9.2 

5.768 

2.41 

89 

11 

11.7 

6.241 

2.23 

113 

12 

12.7 

6.311 

2.62 

102 

13 

13.3 

6.380 

1.90 

116 

14 

13.8 

6.401 

2.91 

127 

15 

14.2 

6.472 

2.85 

108 

16 

14.6 

6.543 

3.62 

111 

17 

14.6 

6.440 

3.04 

94 

18 

12.3 

6.057 

1.61 

105 

19 

8.1 

5.217 

2.06 

119 

20 

6.5 

4.848 

1.64 

120 

21 

5.6 

4.664 

1.79 

111 

22 

4.3 

4.256 

1.39 

104 

23 

2.8 

4.008 

1.80 

103 

24 

2.6 

3.951 

1.81 

98 
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Table  4. 

Meteorological  model 
input  (30-min  average 
data)  for  Davis  13-14 
July  1966  data  set 
(summer). 


Local 

time 

Temperature 

(°C) 

Specific  humidity 
(g/kg) 

Wind  speed 
(m/s) 

Wind  direction 
(°N) 

1 

12.81 

8.08 

1.800 

154 

2 

12.46 

8.03 

1.839 

146 

3 

12.08 

7.96 

1.239 

167 

4 

11.96 

7.92 

0.313 

167 

5 

11.46 

7.69 

0.516 

145 

6 

12.74 

8.14 

1.648 

299 

7 

15.03 

8.48 

1.934 

327 

8 

17.21 

9.12 

1.845 

336 

9 

19.21 

9.63 

1.629 

276 

10 

20.63 

10.11 

1.939 

194 

11 

22.61 

10.22 

1.408 

210 

12 

23.88 

10.11 

2.865 

218 

13 

25.12 

9.92 

3.247 

208 

14 

25.85 

9.33 

5.026 

192 

15 

25.84 

8.67 

5.932 

206 

16 

25.12 

8.24 

5.738 

218 

17 

23.40 

8.16 

5.429 

210 

18 

21.87 

8.10 

5.151 

204 

19 

18.85 

7.72 

4.403 

206 

20 

16.48 

7.30 

4.582 

210 

21 

14.46 

7.47 

3.177 

208 

22 

13.50 

7.46 

2.723 

204 

23 

13.19 

7.50 

3.686 

205 

24 

12.82 

7.49 

3.996 

210 

1 

12.43 

7.67 

3.966 

204 

2 

11.87 

7.47 

3.739 

198 

3 

11.68 

7.44 

3.432 

200 

4 

11.51 

7.35 

3.979 

192 

5 

11.34 

7.27 

2.923 

194 

6 

12.56 

7.6 

1.902 

240 
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4.  Model  Results 


4.1  Surface  Energy  Budget 

Figure  1  shows  computations  for  two  components  of  the  surface  energy 
budget  (net  radiation  and  soil  heat  flux)  for  the  Hay  data  set,  along  with 
field  data.  The  calculated  and  observed  data  appear  well  matched.  As 
these  results  show,  the  model  is  generally  effective  in  computing  radiative 
exchanges  over  near-barren  ground  and  heat  transfer  through  a  dry  soil 
layer. 

On  the  other  hand,  ground  temperature  approximations  (not  shown)  were 
somewhat  underestimated,  as  much  as  1  to  3  °C,  particularly  during  the 
night  and  early  morning  hours  before  sunrise.  Had  the  predicted  surface 
temperatures  been  more  in  line  with  the  observed  data  for  those  times,  the 
calculated  value  for  the  outgoing  soil  flux  would  have  been  higher.  Often, 
discrepancies  of  this  kind  between  data  and  model  output  can  be  mini¬ 
mized  by  a  process  of  iterative  "tuning"  (i.e.,  adjustments  made  within  jus¬ 
tifiable  limits)  of  just  a  few  specific  input  parameters,  such  as  soil  water 
content,  albedo,  or  surface  roughness,  for  example.  In  this  case,  however, 
the  model  was  relatively  unresponsive  to  this  kind  of  tuning. 

Figure  2  shows  estimates  of  the  surface  energy  budget  for  the  Davis  1967 
data  set.  The  model  calculations  agree  remarkably  well  with  the  observed 
data  for  the  net  radiative  flux  and  the  evaporative  flux  above  the  irrigated 
canopy  layer.  The  sensible  heat  and  soil  heat  fluxes,  however,  do  not  agree 
as  well  with  the  observed  data  as  would  have  been  expected  considering 
the  other  computed  results.  In  previous  work  with  this  data  set  (Rachele 
and  Tunick,  1994;  Rachele  et  al,  1996b),  it  was  found  that  increased  wind 
speeds  beginning  from  1300  to  1400  LT  (see  table  4)  were  catalysts  for 
higher  rates  of  evaporation,  drawing  heat  from  the  air  in  the  layer  adjacent 
to  the  canopy,  and  promoting  increased  surface  layer  stability.  In  the  model 
calculations  reported  here,  however,  the  sensible  cooling  of  the  air  did  not 
seem  to  take  place  as  I  would  have  expected.  Figure  2  shows  a  nearly  3-hr 
time  difference  where  =  0.0  W/m^.  The  modeled  excesses  in  the 
sensible  heat  flux  suggests  that  the  calculated  surface  temperatures  (not 
shown)  may  have  been  too  high.  These  surface  layer  processes,  by  the 
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Figure  1.  Model  comparisons  to  observations  of  net  radiative  flux  (left)  and  soil  heat  flux  (right)  for 
Hay  1967  data  set. 
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Soil  heat  flux  (W/m2)  Net  radiative  flux  (W/m2) 


Local  time  Local  time 


Local  time  Local  time 

Figure  2.  Model  comparisons  to  observations  of  net  radiative  flux  (top  left),  sensible  heat  flux  (top 
right),  soil  heat  flux  (bottom  left),  and  evaporative  heat  flux  (bottom  right)  for  Davis  1966  data  set. 

definition  of  microclimate,  affect  the  structure  of  temperature,  wind  speed, 
and  humidity  profiles  through  the  overlying  boundary  layer  (as  I  show  in 
sect.  4.2). 

4.2  Boundary-Layer  Profile  Structure 

From  the  boundary-layer  data  collected  at  Hay,  I  studied  6-  and  18-hr  fore¬ 
casts  in  order  to  evaluate  the  microclimate  model  with  regard  to  mixed 
layer  profile  structure.  Figure  3  shows  modeled  and  recorded  air  tempera¬ 
ture  and  specific  humidity  profiles.  The  model  output  agrees  with  many  of 
the  gross  features  of  the  observed  data:  for  example,  the  simulated  night¬ 
time  temperature  inversions  (top  left)  appear  to  be  captured  by  the  model, 
even  though  the  intensities  of  the  near-surface  and  upper-level  gradients 
differ.  The  computed  temperature  profiles  for  the  daytime  case  (top  right) 
also  appear  to  be  in  agreement  with  observed  data;  however,  the  tempera- 
tiure  profile  gradients  at  the  top  of  the  boundary  layer  are  clearly  offset  (see 
fig.  4).  On  close  inspection,  the  model  and  observed  daytime  temperature 
surface  gradients  are  offset  as  well. 


*As  a  matter  of  convention,  for  the  net  radiative  and  soil  heat  flux,  positive  values  represent  energy  transferred  down¬ 
ward,  while  negative  values  represent  energy  transferred  away  from  the  soil-plant-atmosphere  interface.  The  opposite 
is  an  accepted  convention  for  the  sensible  and  evaporative  heat  fluxes. 
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Figure  3.  Model  comparisons  to  observations  of  profiles  of  boundary-layer  temperature  (top)  and 
specific  humidity  (bottom)  for  Hay  1967  data  set. 


Figure  4.  Model 
comparison  to 
observations  of  depth 
of  planetary 
boundary  layer  for 
Hay  1967  data  set. 


The  simulated  specific  humidity  profiles  in  figure  3  (bottom)  show  offsets 
in  water  vapor  content  through  the  upper  and  lower  levels  of  the  bound¬ 
ary  layer,  in  contrast  to  the  observed  data.  The  observed  data  indicate  a 
characteristic  drying  out  of  the  lower  atmosphere  due  to  vertical  mixing. 
Alternatively,  the  modeled  profiles  indicate  tendencies  toward  an  overall 
increase  of  boundary  layer  moisture  over  this  period  of  time. 
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By  definition,  mixing  depths  are  determined  as  a  function  of  boundary- 
layer  growth  and  decay.  These  processes,  in  general,  are  driven  by  the 
radiative  exchanges  that  govern  the  surface  heat  and  moisture  budget. 
Therefore,  discrepancies  between  predicted  and  observed  inversion 
heights  are  most  often  due  to  the  inaccuracy  of  the  predicted  surface-layer 
gradients.  Otherwise  influences  on  boundary-layer  growth  must  be  attrib¬ 
uted  to  external  and  upper-level  forcings,  such  as  advection  or  changes  in 
the  large-scale  horizontal  pressure  gradient.  The  unmatched  inversions 
shown  in  figure  3  are  then  most  likely  the  result  of  overpredicted  surface- 
layer  heating.  This  was  found  to  be  the  case  in  an  earlier  study  that  was 
repeated  by  Pielke  and  Mahrer  (1975)  using  the  prediction  expression  of 
Deardorff  (1974)  (eq  (4),  sect.  2,  this  report).  Values  of  calculated  mixed- 
layer  heights,  shown  in  figure  4,  compared  to  those  determined  from  ob¬ 
servations  are  given  in  table  5. 

Figure  5  shows  model-derived  wind  speed  and  wind  direction.  The  ID 
model  predicts  an  observed,  nighttime  wind  speed  maximum  (left)  very 
poorly.  It  also  fails  to  match  fluctuations  in  the  observed  daytime  wind 
speed  profile  (right)  although  on  average,  the  mean  differences  are  less 
than  1  m/s. 

Diurnal  accelerations  and  decelerations  in  wind  speed  are  normally  due  to 
changes  in  mesoscale  pressure  gradients  (over  a  range  of  20  to  200  km). 
These  gradient  changes  are  brought  about  by  changes  in  temperature  near 
the  ground,  which  increase  or  decrease  aerodynamic  drag  and  stress 
(McNider  and  Pielke,  1981).  To  model  these  processes  adequately,  one 
would  need  time-dependent,  3D,  horizontal  and  vertical  transfer  equa¬ 
tions  for  momentum,  mass,  and  temperature  for  each  point  in  a  user- 
defined  grid.  From  the  results  of  this  study,  it  becomes  apparent  that  the 
ID  approach  is  not  sufficient. 

The  forecasts  of  wind  direction  profile  structure  shown  in  figure  5  (bottom) 
are  also  poor  through  the  forecast  period,  even  though  mean  differences 
may  appear  to  be  on  the  order  of  ±15°.  It  is  known  that  wind  directions  are 
probably  the  most  difficult  variable  for  meteorological  simulation  models 


Table  5.  Calculated 
heights  of  planetary 
boundary  layer 
compared  to 
observations. 


Local  time 

Height  of  bovindary  layer  (m) 

Observed  data^ 

Model 

0800 

— 

100.00 

0900 

110.00 

107.36 

1000 

— 

175.20 

1100 

— 

439.79 

1200 

1000.00 

1018.05 

1300 

— 

1373.03 

1400 

— 

1448.27 

1500 

1160.00 

1464.05 

1600 

— 

1468.23 

1700 

— 

100.00 

1800 

1240.00 

100.00 

‘^Estimated  from  observed  potential  temperature 
and  specific  humidity  profile  data  (Pielke  and 
Mahrer,  1975). 
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to  predict.  Events  on  a  number  of  different  scales,  from  frontal  passages  to 
local  urban  or  rural  land  use,  affect  wind  direction.  Assessments  of  pre¬ 
dicted  wind  direction  have  shown  tremendous  scatter  compared  to  obser¬ 
vations  (see  fig.  13.2  in  Gross,  1994).  However,  next  to  temperature  and 
changes  in  temperature  with  height,  wind  speed  and  direction  are  signifi¬ 
cantly  important  in  assessing  the  outdoor  acoustic  field  that  results  as  a  de¬ 
rivative  of  a  microclimate  forecast.  The  last  sections  of  this  report  discuss 
the  assessment  of  outdoor  acoustic  fields  in  greater  detail. 


Figure  5.  Model  comparisons  to  observations  of  profiles  of  boundary-layer  wind  speed  (top)  and  wind 
direction  (bottom)  for  Hay  1967  data  set. 
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5.  Application  to  Atmospheric  Acoustics 

5.1  Calculated  Boundary-Layer  Sound-Speed  Profiles 

First-order  approximations  of  sound  speed  in  the  absence  of  wind  can  be 
derived  from  the  expression  in  Pierce  (1981)  given  as 


where  y  is  the  ratio  of  the  specific  heats  of  air  for  constant  pressure  and 
constant  volume  modified  to  account  for  water  vapor  content  so  that  y=  (7 
+  h)/ (5  +  h),  where  h  is  the  fraction  of  water  molecules  in  the  air,  R  is  the 
universal  gas  constant  (8314.16  kg"'  K“^),  T  is  air  temperature  in  kelvins, 
and  M  is  molecular  weight,  calculated  as  M  =  29  -  llh.  However,  the  varia¬ 
tions  of  sound  speed  caused  by  even  extreme  changes  in  humidity  are 
minimal  (approximately  1  to  2  m/s)  and  can  generally  be  ignored  (Noble 
and  Marlin,  1995).  The  effect  of  winds  on  sound  speed  is  expressed  by  a 
vector  relation.  Equation  (6)  can  be  rewritten  as 

C,ff (z)  =  CiT)  +  U cos  Br)  .  (7) 

where  C(T)  is  the  sound  speed  in  the  absence  of  wind  at  temperature  T,  U 
is  the  mean  horizontal  wind  speed  (in  meters  per  second),  6^,  is  the  wind 
direction  (conventionally  interpreted  as  the  direction  from  which  the  wind 
originates),  6^^-  k  is  interpreted  as  the  along-wind  direction,  and  Br  is  the 
azimuth  or  bearing  of  an  outlying  receiver  from  a  given  acoustic  source. 
The  environmental  effects  on  the  propagation  of  acoustic  energy,  as  re¬ 
viewed  by  Noble  (1992),  are  related  to  the  height-dependent  variations  of 
the  sound-speed  profile,  which,  in  general,  cause  (1)  upward  refracting 
conditions,  when  the  sound-speed  gradients  through  the  first  0  to  400  m  of 
the  mixed  layer  are  negative  (normally  occurring  during  the  day),  or 
(2)  downward  refracting  conditions,  when  the  sound-speed  gradients  are 
positive  (normally  occurring  at  night). 

Figure  6  shows  the  derived  boundary-layer  sound-speed  profiles  for  the 
nighttime  and  daytime  cases.  There  appears  to  be  reasonably  good  agree¬ 
ment  between  the  sound-speed  profiles  derived  from  model  calculations 
and  those  derived  from  the  field  data,  particularly  for  the  daytime  case 
shown.  Differences  between  the  results  are  principally  related  to  the  com¬ 
bined  accuracy  of  simulated  gradients  in  temperature  and  wind-speed 
profile  structure.  (It  is  only  by  a  matter  of  convention  that  the  calculations 
in  the  bearing  Bsu,face  (left)  are  called  upwind,  while  those  made  in  the  bear¬ 
ing  Bsu^ce  - 180°  (right)  are  called  downwind.) 
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Figure  6.  Derived  boundary-layer  sound-speed  profiles  for  nighttime  (top)  and  daytime  (bottom) 
atmospheric  conditions. 


5.2  Approximations  of  Short-Range  Outdoor  Acoustic 
Attenuation 

Meteorological  profiles  from  observed  Hay  field  data  and  those  forecast  by 
the  model  in  this  study  were  used  as  input  to  the  acoustic  propagation  nu¬ 
merical  code  WSCAFFIP  (Windows  (version)  Scanning  Fast  Field  Pro¬ 
gram),  which  assesses  environmental  effects  on  short-range  acoustic 
attenuation.  WSCAFFIP  provides  attenuation  levels  with  range  and  azi¬ 
muth  for  a  given  geometry  and  frequency  (Noble,  1996;  Noble  and  Marlin, 
1995).  The  acoustic  propagation  algorithms  in  WSCAFFIP  attempt  to  rep¬ 
resent  the  effects  of  atmospheric  refraction,  diffraction,  absorption,  and 
reflection  or  ground  impedance.  Figure  7  shows  the  WSCAFFIP  results  for 
the  upwind  (top)  and  downwind  (bottom)  nighttime  cases.  Ahead  of  the 
first  segment  of  the  data  (i.e.,  x  <  650  m),  the  calculated  attenuation  levels 
appear  to  be  well  matched.  Further  out  in  range,  differences  in  the  derived 
sound-speed  profile  gradients  have  resulted  in  either  sound-limited  areas 
(shadow  regions)  or  areas  of  refocusing  due  to  refraction  (i.e.,  the  local 
minima  and  maxima).  As  a  rule,  the  more  positive  the  sound-speed  gradi¬ 
ent  (see  fig.  6  (top)),  the  greater  the  effects  of  a  downward  refracting 
atmosphere,  especially  in  the  case  of  calculations  made  in  the  downward 
direction. 
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Figure  7.  Short-range 
acoustic  attenuation 
upwind  (top)  and 
downwind  (bottom) 
predicted  by 
WSCAFFIP 
numerical  code  for 
nighttime 
atmospheric 
conditions. 


In  contrast,  the  attenuation  levels  for  the  daytime  case  (fig.  8)  show  the  ef¬ 
fects  of  upward  refracting  atmospheric  conditions.  The  amount  of  attenua¬ 
tion  appears  to  be,  overall,  greater  than  in  the  nighttime  case.  Soimd  levels 
are  shown  to  be  in  good  agreement  to  a  distance  of  approximately  1750  m 
in  range.  Beyond  1750  m,  offsets  in  the  profiles  predictions  have  affected 
the  degree  of  separation  in  the  data,  as  seen  in  both  figures  in  the  range 
from  2000  to  4000  m. 
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Figure  8.  Short-range 
acoustic  attenuation 
upwind  (top)  and 
downwind  (bottom) 
predicted  by 
WSCAFFIP 
numerical  code  for 
daytime  atmospheric 
conditions. 
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6.  Conclusions 

In  the  work  reported  here,  I  used  micrometeorological  field  data  in  an 
effort  to  evaluate  an  atmospheric  microclimate  model.  I  set  out  to  deter¬ 
mine  if  it  was  possible  to  forecast  ID  profiles  of  the  atmospheric  boundary 
layer  sufficiently  well  that  the  corresponding  outdoor  acoustic  field  could 
be  reliably  calculated. 

Specifically,  I  studied  a  soil— plant— atmosphere  model  using  two  dissimilar 
sets  of  experimental  data;  Hay  1967  and  Davis  1966.  Model  calculations  of 
the  surface  energy  budget  were  compared  to  field  observations.  Derived 
meteorological  profiles  for  the  overlying  boundary  layer  were  also  com¬ 
pared  to  data.  Using  these  modeled  results,  I  derived  outdoor  sound-speed 
profiles.  I  also  used  the  modeled  results  as  input  to  the  short-range  acous¬ 
tic  propagation  numerical  code  WSCAFFIP. 

Some  general  conclusions  with  regard  to  the  evaluation  of  this  atmo¬ 
spheric  microclimate  model  follow. 

1.  Quantitative  assessments  of  the  surface  energy  budget  over  land  princi¬ 
pally  rely  upon  equations  sensitive  to  variations  in  soil  wetness  and  sur¬ 
face  temperature,  since  these  are  probably  the  two  factors  most  sigmficant 
determining  the  partitioning  of  the  soil,  evaporative,  and  sensible  heat 
fluxes. 

2.  As  shown  by  the  results  of  the  surface  energy  budget  calculations  with  the 
Hay  1967  field  data  used  as  input,  the  model  is,  overall,  effective  in  com¬ 
puting  radiative  exchanges  over  sparsely  covered  ground  and  heat  trans¬ 
fer  through  a  dry  soil  layer. 

3.  Discrepancies  between  data  and  model  output  can  normally  be  minimized 
by  a  process  of  iterative  "tuning"  of  just  a  few  specific  input  parameters, 
such  as  soil  water  content,  albedo,  or  surface  roughness.  For  the  case  with 
the  Davis  1966  data  used  as  input,  however,  the  model  was  relatively  unre¬ 
sponsive  to  this  kind  of  tuning.  Therefore,  although  the  radiative  and 
evaporation  flux  calculations  for  the  irrigated  canopy  were  in  good  agree¬ 
ment  with  the  field  data,  the  sensible  and  soil  heat  fluxes  agreed  less  well 
with  the  observed  data  than  I  had  expected. 

4.  The  ID  forecasts  of  meteorological  profile  structure  in  the  atmospheric 
boundary  layer  were  generally  very  poor,  especially  with  regard  to 
windspeed.  This  level  of  model  performance  was  disappointing.  The  ID 
approach  appears  to  be  insufficient  to  effectively  represent  effects  of  sur¬ 
face  roughness,  wind  shear,  and  the  exchanges  of  heat  and  moisture  from 
one  level  to  the  next. 

5.  Quantitative  assessments  of  the  outdoor  acoustic  field  rely  upon  accurate 
representations  of  the  local  meteorology,  since  variations  in  the  profile  gra¬ 
dients  of  soimd  speed  are  determined  from  these  data.  Therefore,  the  dis¬ 
crepancies  in  the  forecasts  of  short-range  acoustic  attenuation  with  range 
presented  in  this  study  were  most  likely  due  to  sufficiently  significant  off¬ 
sets  in  the  profiles  of  sound  speed,  particularly  in  the  0  to  400  m  layer 
above  ground. 
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7.  Recommendations 


The  components  of  the  atmospheric  surface  energy  budget,  by  definition, 
reflect  the  rates  of  heating  and  cooling  near  the  ground.  The  transfer  of 
incident,  reflected,  and  emitted  energy  at  the  surface  can  be  approximated 
as  functions  of  cloud  cover,  terrain  roughness  characteristics,  surface  al¬ 
bedo,  canopy  leaf  area,  soil  water,  and  soil  thermal  capacities.  As  an  inte¬ 
gral  part  of  the  microclimate,  the  surface  energy  budget  acts  as  the  princi¬ 
pal  influence  on  the  development  of  pressure,  temperature,  and  humidity 
gradients  on  all  scales.  As  a  result,  most  of  the  turbulence  and  diffusion  of 
heat,  momentum,  and  moisture  in  the  boundary  layer  can  be  attributed  to 
surface  layer  energy  exchanges.  They  are  inherently  multidimensional 
processes. 

In  recent  years,  however,  several  efforts  have  relied  on  ID  meteorology  for 
defense-related  applications,  in  particular,  outdoor  acoustics.  Both  models 
and  measured  field  data  have  been  manipulated  in  efforts  to  create  reliable 
product  codes.  Unfortunately,  at  the  conclusion  of  each  initiative,  it  is 
found  that  better  assessments  or  forecasts  of  wind,  temperature,  and  hu¬ 
midity  fields  require  further  research  (into  models  that  can  account  for 
changes  in  the  state  variables  over  time,  both  horizontally  and  vertically, 
from  within  as  well  as  from  outside  the  modeling  area  of  interest).  The  ID 
approach  applied  to  areas  10  to  80  km  across  has  been  insufficient,  espe¬ 
cially  with  regard  to  estimates  of  wind  speed  and  wind  direction.  In  appli¬ 
cation,  only  limited  successes  have  been  realized.  Therefore,  because  of 
their  explicit  dependence  on  temperature  and  wind  profile  structure, 
acoustic  attenuation  models,  especially,  need  to  incorporate  3D  meteoro¬ 
logical  data  (or  forecast  fields).  Their  reliability  in  the  battlefield  environ¬ 
ment  would  then,  in  my  opinion,  greatly  improve. 
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